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Abstract 

The asymptotic time derivative of the total dipole signal is proposed as an useful observable 
to investigate on Isospin equilibration phenomenon in multi-fragmentation processes. The study 
has been developed to describe charge/mass equilibration processes involving the gas and liquid 
"phases" of the total system formed during the early stage of a collision. General properties of this 
observable and the links with others isospin dependent phenomena are discussed. In particular, 
the 40 C7 + 28 Si system at 40 MeV/nucleon is investigated by means of semiclassical microscopic 
many-body calculations based on the CoMD-II model. The study of the dynamical many-body 
correlations produced by the model also shows how the proposed observable is rather sensitive to 
different parameterizations of the isospin dependent interaction. 
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I. INTRODUCTION 

An interesting subject related to Heavy Ions Isospin physics [l| is the process leading 
to the equilibration of the charge/mass ratio between the main partners of the reaction 



as well described in Ref. 



I]- 



The so-called "isospin diffusion" phenomenon is the re 



evant 



mechanism acting between the reaction partners in binary processes 
particular, in the collision of the 124 and 112 Tin isotopes at 50 MeV/nucleon 3j, evidence of 
partial equilibrium in the charge/mass ratios of the quasi-projectile and quasi-target has been 
deduced through the study of the iso-scaling parameters related to the isotopic distributions. 
In this case dynamical calculations, based on the Boltzmann-Uehling-Uhlenbeck model s), B] , 
show that the degree of equilibration depends on the behavior of the symmetry potential 
U T as a function of the density. 

In this work we want to extend the study of the isospin equilibration processes, looking 
at the system in a global way, by using the following quantity: V(t) = J2f=i At a 
microscopic level, the sum on the index i is performed on all the Z to t protons (bound and 
free) of the system. V (t) corresponds, apart from the elementary charge e, to the time 
derivative of the total dipole of the system. The velocities ~v\ are computed in the center of 
mass (cm.) reference frame. We note that, as due to the total momentum conservation, the 
global effect related to the motion of the neutral particles (bound and free) is also implicitly 
contained in V^(t). Therefore we can expect, in a quite general way, a peculiar dependence 
of the behavior of V^(i) from the iso- vectorial interaction. 



Several studies were based on this dynamical varia 



jle to^ describe pre-equilibrium Giant 
and references therein ). Var- 



Dipole Resonance (GDR) 7-ray emission (see Ref. Q, flQ 
ious reasons suggest us to use the same variable to also describe isospin equilibration in 
complex processes. The basic guide line starts just from these studies on pre-equilibrium 
radiative emission in essentially binary processes induced by charge/mass asymmetric part- 
ners. The ensemble average of the time derivative of the total dipole, in the initial and final 



stages of the collision, can be approximated through the so-called Molecular component [131 ]: 



(?) « (?m) = \{»ptWt) - (Pp))(v^) (1) 
where \ip? is the reduced mass of the projectile or quasi-projectile (P) and target or quasi- 
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target (T) of the binary system, {3 P>T = Np '^~^ P ' T represents the associated relative neutron 
excess (Np t T, Zp t T, Ap t T are the neutron, proton and mass numbers respectively), vpt in- 
dicates the relative velocity During the interaction stage of the collision, if the partners 
have exchanged charge and mass for enough time in such a way to equilibrate, on average, 
the charge/mass ratio, then the final value of (V) will be zero. In this case the yield of 
the pre-equilibrium dipolar 7-ray emission, which satisfies the well known selection rule on 
the isospin quantum number T (AT = 1, no isospin mixing) of the involved intermediate 
states, reaches the maximum value. For more complex mechanisms, as the ones leading to 
the substantial stopping of the two partners and to the formation of one hot source, the 
dipolar signal related to the other particles and fragments can not be neglected. However, 
also in these cases we can use the same definition of isospin equilibration through the con- 
dition (V) = 0. As we will show in the following, this definition leads to other conditions 
concerning the production of differential flow between neutral and charged nucleons 
and/or the relative neutron excesses of the main produced fragments and free particles. 

In the following we briefly discuss some simple examples and properties of ( V ) = aiming 
to illustrate the information which is potentially contained in this quantity and the relation 
with charge/mass equilibration processes: 

- (i) By taking into account only the effects associated to the strong interaction, af- 
ter the pre-equilibrium stage, starting from the time t pre , when a second stage charac- 
terized by an average isotropic or symmetric emission of the secondary sources (statisti- 
cal equilibrium) takes place, the ensemble average of V (t) satisfies the following relation: 



(V(t pre )) = (V(t > t pre )) = ( V*) U|- The average value of this dynamical variable at t pre is 
invariant with respect to statistical processes and therefore the value of (v) is determined 
only by the complex dynamics which characterizes the early stage of the collision, when fast 
changes of the average nuclear density are expected. In particular, (V ) can be expressed 
as a function of the charge Z, mass A, average multiplicity (mz,A) and the average value of 
the mean momentum (P z,a) per unit of mass (expressed in fmT 1 ) of the detected particles 
in the generic event: 

^)=Y.-Mz,a)^za)C Z 4 (2) 

Z,A A F 

qZ,a _ {mz,A~P z,a) 
~P " (P ZtA )(m ZtA ) 



C-^ is the correlation function between the multiplicity and the mean momentum. This 
correlation function plays a key role for the invariance property and therefore requires for 
an event by event analysis in which many-body correlations can not be neglected. For 
symmetry reasons, (V) lies on the reaction plane. It is directly linked with a weighted 
mean of the charge/mass ratio, as Eq.(3) suggests. It also takes into account the average 
isospin flow direction through the momenta ~P z A- The long range Coulomb interaction 



can produce differences between the value of (V(t pre )} and the observed asymptotic value. 
These changes, however, are rather small at the involved energies (see Sec. Ill) and can be 
evaluated with the necessary precision by taking into account the corrections due to the 
Coulomb repulsion up to the asymptotic stage. 

-(h) In the general case, we find attractive the following decomposition: (v) = ( V a) + 
(V l) + (vgl) where (v g) and (v l) are the average dipolar signals associated to the gas 
"phase" (light charged particles) and to the "liquid" part 16] corresponding to the motion of 
the produced heavy fragments. The signal ( V*gl) is instead associated to the relative motion 
of the two "phases". By supposing, for simplicity, the gas "phase" formed by neutrons and 
protons, (v) can be further decomposed as: 

(?) = ( MLiM ^ + ( ^(gpM ^^) + { v r , L ) ( 4) 

In the above expression the first term represents the contribution related to the proton- 
neutron relative motion of the gas "phase" composed of Aq nucleons, expressed through the 
relative velocity lf^ N ; the second term concerns the relative velocity lf C m,LG between the 
centers of mass of the "liquid" complex and the "gas", [Iql is the reduced mass associate 
to the two sub-systems; the last term represents the contribution produced by the relative 
motion of the fragments. A similar expression can be obtained including other light particles 
in the gas "phase" as , for example, light Intermediate Mass Fragments (IMF). From this 
decomposition we can see how the isospin equilibration condition ((V) = 0), for the total 
system, requires a delicate balance which depends on the average neutron excess of the 
produced "liquid drops" (Pl), on the one associated to the gas "phase" and on the 

relative velocities between the different parts. To enlighten the role played by some of the 
terms reported in Eq.(4), we can discuss the idealized decay of a charge/mass asymmetric 
source through neutrons and protons emission (or the case in which the liquid drops are 
produced through a statistical mechanism (V r ,L) — 0)- This example can also schematically 



describe the case of the complete stopping in heavy ion collisions. Moreover, for simplicity, we 
can consider uncorrelated fluctuations between the velocities, masses and neutron excesses. 
In absence of pre-equilibrium emission or for identical colliding nuclei the second term of 
Eq.(4) is zero ((if ' C m,LG) = 0), and the isospin equilibration requires a neutron or proton 
gas "phase" or absence of relative neutron- proton motion. For non-identical colliding nuclei, 
if pre-equilibrium emission exists, then {l?cm,LG) 0- m this case, if (Pa) (Pl), due, 
for example, to the isospin "distillation" phenomenon , the first term has to be necessarily 
different from zero and it will contribute to the neutron-proton differential flow (see also 
Sect. III). We remark that both the isospin "distillation" phenomenon 13] and the production 
of neutral-charged particles differential flow [15j depend in a peculiar way from the iso- 
vectorial forces. 

Therefore, according to our description, the understanding of the isospin equilibration pro- 
cess for the total system requires the gas "phase" contribution to be taken into account. 
This term can be regarded as a kind of " dissipation" with respect to the system formed by 
the liquid part. In this work, as an example, we will discuss the results obtained through the 
Constrained Molecular Dynamics-II approach (CoMD-II) 18|,ll9|] applied to the charge/mass 
asymmetric A0 Cl+ 28 Si system at 40 MeV/nucleon. The study is performed by using different 
options for the iso-vectorial potential term. 

Before to show the results of our calculations, in the following section we briefly recall the 
way in which the isospin dependence of the nuclear interaction is introduced in the CoMD-II 
model. 



II. SYMMETRY INTERACTION AND CORRELATIONS 

According to the results shown in Ref . [rsl. I20I . starting from a Skyrme type two-body mi- 
croscopic interaction, the two-body effective potential in CoMD-II model can be expressed 
through the nucleon-nucleon overlap integral p lJ = J / ^rjd^r^^ — r^)p l fp . ~fj, f\ repre- 
sent the nucleon spatial coordinates. p l is the Gaussian distribution in the coordinate space 
related to the generic nucleonic wave-packet. The microscopic iso-vectorial interaction for 
the Stiff2 option is given by the following expression: V T = ^ J2^tj=i(^S TitTj — l)5(~r*j — ~r*j), 
A is the total mass number, Tj represent the generic third nucleonic isospin component and 
Po is the one-body density at the saturation point. The coefficient a = 72MeV deter- 
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mines the strength of the iso-vectorial interaction at the saturation density. As shown in 
Ref. [20] the structure of V r can be obtained by taking into account that the two-body nu- 
clear forces, in S wave, around the ground state (g.s.) density, are less attractive in isospin 
triplet states (T = 1) with respect the singlet states (T = 0). From the above expressions it 
results that the associated effective interaction U T (after the convolution with the nucleonic 
wave-packets) can be expressed as a function of the average overlap integrals per couple of 
neutrons (nn), p nn , protons (pp) p pp , and neutron-proton (np) p np . As discussed in Ref.j^j 
for small asymmetries we can assume p nn = fP p = - — 1^— = p. To characterize the dif- 
ferences associated to the nucleon-nucleon dynamics, at a two-body level, we can introduce 
the correlation coefficient a in such a way p np = (1 + a) p. a depends on both p and the 
asymmetry parameter (3. Results on nuclear matter simulations 20] show that the behavior 
of a as a function of ft can be approximated for moderate asymmetries by a parabolic law. 



In this case, eqs.(6,7) of Ref. 20j give the following expression for the effective iso - vectorial 



potential in the Non Local (N.L.) approximation U T N L 

U t n .l. = ^pAWisW + \a - a')(3 2 - \a \ (5) 

a = iQ^y-o (6) 

(3 A terms are neglected in the previous expression. a = a (f>, ft — 0) represents the correlation 
coefficient related to the difference in the dynamics of the np couples with respect to the nn 
and pp ones for symmetric nuclear matter. 

It depends on the average overlap integral per couple of nucleons p which reflects the 
degree of compression, s = ^ J2i^j P %)i is associated to the total overlap integral per nucleon. 
F' is a form factor which modulates the changes of the iso-vectorial interaction as a function 
of the average overlap integral s. For the Stiff 1 option we use F' = - 2s +s , for the Stiff2 
case F' = 1 and for the Soft option F' = (^i.) 1 / 2 . These form factors correspond to the 

nn 

following values of the K sym parameter as defined in Ref. [14|, [15( : 94 MeV for the Stiff 1 
option, -27 MeV for the Stiff2 option and -88 MeV for the Soft one. This choice corresponds 
also to an Sq value [14j of about 40 MeV . From Eq.(5) we note that in our approach the 
iso-vectorial forces generates, beyond the (3 dependent potential, also another iso-vectorial 
density dependent term, independent on (3 and proportional to the degree of correlation a 



evaluated for symmetric systems. As discussed in 



2d, 



21], at small asymmetries, this term 



determines the high sensitivity of the experimental observables to the different functional 
forms of F' . The finite value of a, which is of the order of 15%, apart from the Pauli principle 
and the Coulomb interaction, is strongly affected by the iso-vectorial interaction itself. In 
fact the neutron-proton couples, at variance with the other ones, suffer the more attractive 
singlet interaction. 

For clarity, we have to note that the np correlations generated by similar kind of forces 
represent a rather interesting subject for the description of the extra-bind energies in iV ~ Z 
nuclei. These studies are performed with static sophisticated quantal approaches. The np 



pairing correlations, for exampl e, p 



energy in symmetric nuclei 22 



ay a main role to microscopically describe the Wigner 



241 ] . However, apart from effects induced by the Pauli 



principle in phase-space (as treated in our approach) and from the usage of an iso-vectorial 
potential which reflects a two-body interaction dependent on isospin quantum numbers, the 
description and the propagation of the many-body correlations in our model calculations 
are essentially classical. The np correlations discussed in the present work can represent 
a dynamical classical analogous which includes Pauli principle and isospin effects of more 
complex quantal effects. Therefore the effects here discussed can not be compared with the 
ones able to go beyond this classical analogous and having therefore a pure quantal nature. 
In fact our approach can not describe effects that, for example, could specifically arise from 
details of the single particles wave functions or also from details concerning the propagation 
and symmetries of the Slater determinant structure of the many-body wave functions, which 



25 



26]. 



instead characterizes the Fermionic molecular dynamics approaches 
The case corresponding to vanishing values of the correlation a represents, in our frame- 
work, the so-called Iso-vectorial Mean Field Approximation (I.M.F.A.). In this case the 
average overlap integrals per couple of nucleons related to neutron-neutron, proton-proton 
and neutron-proton interactions have the same values (p nn = fp p = p np ) and the iso-vectorial 
interactions generate only the usual symmetry potential term which depends on {3 2 . 
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III. CALCULATION RESULTS 



A. An example: the 40 C7 + 28 Si system at 40 MeV/nucleon 

Now we discuss, as an example, the results concerning the isospin equilibration process 
for the 40 CI + 28 Si system at 40 MeV/nucleon. In Fig. 1 we show the average total dipolar 
signals evaluated through CoMD-II calculations along the z beam direction (V z ) and along 
the impact parameter direction x, (V x ), respectively. The reference frame is the cm. one. 
The impact parameter b is equal to 3 fm, in panels (a) and (c) and 1.5 fm in panels (b) 
and (d). In Fig. 1(a) and Fig. 1(c) the average dipolar signals are shown for the first 
140 fm/c. Different lines refer to different iso-vectorial potentials, according to Ref. 18]. 



In the first 150 fm/c in all the cases wide oscillations exist. They are responsible for the 



pre-equilibrium 7-ray emission 
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11| . The damped oscillations converge towards smaller 
and almost constant values. This can be seen in Fig. 1(b) and Fig. 1(d) in which the 
dynamical evolution is followed from 110 fm/c up to 300 fm/c. The inclusion of corrections 
due to Coulomb interaction at longer time is of the order of some percent. The time interval 
in which the almost stationary behavior is reached is related to the lifetime of the coherent 
dipolar collective mode and it is strictly linked with the average time for the formation of 
the main fragments and pre-equilibrium emission. Fig. 2 shows the asymptotic values of 
the z and x components of the average dipolar signal, evaluated for different reduced impact 
parameters b r = (b max ~ 7.5fm) and different interaction options. The figure enlightens 
the sensitivity of the (v) observable to the iso-vectorial interaction. 

In the following we would like to discuss the behaviors of the different components of the 
total dipolar signal along the z direction as a function of the impact parameter according to 

the scheme depicted in the previous section and related to Eqs. (1,4). In Fig. 3(a) the ratio 

(v z ) 

Rgl = |y§y between the z asymptotic components of the dipolar signal associated to the light 
particles (V<4) and the one corresponding to the two biggest fragments (V£) is plotted as a 
function of the reduced impact parameter b r for the Stiff2 option. The other options produce 
similar results. For peripheral collisions we see the predominance of the signal carried by 
the two biggest fragments (\Rgl\ < 1)- In these cases the liquid part is dominated by 
the so-called "molecular" component Vm (see Eq.(l)). It has a negative sign because the 
quasi-projectile has the largest neutron excess with respect to the quasi-target. The "gas" 
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component produces a small negative value which indicates, on average, a dominance of 
the pre-equilibrium emission of nucleons from the target light partner. For smaller impact 
parameters the molecular component reduces its value because of the increasing stopping 
(reduction of vpr in Eq.(l)) and as a consequence of the charge/mass equilibration between 
the two partners. The asymptotic behavior of Rgl between b r = 0.6 — 0.75 is due to the 
change of sign of the "liquid" component which evolves in a rather continuous way passing 
through zero. The "liquid" part has a positive sign for more central collisions, because now 
the process can be roughly described as one charged source emitting pre-equilibrium charged 
particles mostly on the target side. The average recoil of the source in the projectile direction 
produce the positive value of the dipolar signal. The absolute value of Rgl increases for 
more central collisions due to the increasing contribution of the pre-equilibrium particles. 
As discussed in the introductory section, when the collision partners are not identical nuclei 
the isospin equilibration process produces a neutron-proton differential flow contribution if 
the neutron and proton "gases" have different cm. velocities. These pre-conditions are 
verified for the studied collision. 

We briefly recall that the neutron-proton differential flow is an observable that has been 
proposed to investigate symmetry potential effects [l4, [lsj]. It is expressed as the average 
difference between the transverse momenta of the emitted free neutrons and protons having 
rapidity y'\ 

1 A{y ' ] 

Fn P = w iVix{y') (7) 

with Wi = 1 for neutrons and Wi = —1 for protons. For b = 3 f m and for the Stiff 1 and 
Soft options, in Fig. 3(b) we show the neutron-proton differential flow F np , expressed in c 
units, as a function of the particles rapidity, y, normalized to the projectile one y' beam - The 
rapidity values are evaluated in the cm. of the total system reference frame. The dashed 
vertical lines indicate the projectile and target reduced rapidity. 

From the figure we can understand that, by averaging on the rapidity, the neutron-proton 
transversal velocity has a negative value. This reflects the average "bending" of the relative 
motion between the cm. of the emitted neutrons and protons, through the half-plane 
opposite to the impact parameter direction. In Fig. 3(c) we display the average number of 
nucleons as a function of the rapidity for the two options. In both the cases the collision 
produces rather similar "stopping". The results shown in Fig. 3(b) can be compared with 



9 



the calculations displayed in Fig. 3(d) obtained by subtracting, event by event, the cm. 
relative neutron-proton motion related to the "gas" phase. As can be seen, similarly to the 
case of identical nuclei, this correction restores (within the errors associated to the statistics 
of simulations) the almost specular behavior of F np with respect to the rapidity axes. The 
correction acts also along the beam direction. The largest absolute values of F np around 
\y\ ~ 1 for the Soft case is mostly due to the largest neutron excess of the "gas" phase with 
respect to the one obtained for the Stiff 1 option. Finally, in Fig. 3(e) we show the vectors 
associated to the average proton-neutron relative motion if™ f° r the two options. These 
quantities determine the correction on the differential flow due to the isospin equilibration 
process and are associated to the "gas" component of the total dipolar signal (see first term 
of Eq.(4)). 

B. Sensitivity to the different options of the iso-vectorial interaction 

In the following we want to discuss in some detail the sensitivity of the dipolar signal to 
different options concerning the iso-vectorial interaction including also the role played by 
the finite value of the correlation coefficient a. 

According to the nature of the discussed quantity, we can expect a particular sensitivity 
to the iso-vectorial forces because it contains implicitly information on the global relative 
motion of the charge particles with respect the neutral ones. In Fig. 4(a) we show as 
a function of the reduced impact parameter b r the asymptotic values of (V x ) and (V z ). 
The arrow indicates the direction of increasing impact parameters corresponding to the 
marked points with a step Ab r = 0.1. Different symbols represent different options. In 
the region of the bending of the lines, around b = 3 fm, we observe the greater sensitivity 
to the iso-vectorial interactions. This result is particularly evident by studying the ratios 
R = (V X )/(V Z ). In Fig. 4(b) we in fact show the relative change r = ^ between couples 
of different options. We can see that for b r less than about 0.6 large changes are predicted 
according to the different shapes of the form factor F' . This impact parameter region 
produces a substantial stopping of the incident nuclei and is clearly dominated by large 
overlap between projectile and target which gives rise to processes changing from incomplete 
fusion reactions to IMF production . The region of intermediate impact parameters show 
the higher sensitivity when the mechanism evolves with respect to the essentially binary 
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processes which take place at the higher impact parameters. For b r greater than 0.6, in fact, 
the sensitivity is strongly reduced. 

In particular, according to what previously observed (see for example point (ii) of Sec. I), we 
have evaluated the partial contributions (V£) and (V£) related to the two main fragments. 
As an example, for b = 3 fm and for the Stiff2 option, the "liquid" asymptotic values 
are (V£) = -0.120c and (V£) = 0.162c while the total contributions are (V x ) = 0.044c 
and (V z ) = —0.027c. Therefore, it results that the contributions carried by the two main 
fragments only partially contribute to the isospin equilibration process. The remaining part 
("gas"), which in this case we have associated to particles and to the IMF, generates a term 
with opposite sign and similar strength for both directions. It contributes in a decisive way 
to the global equilibration process. For the same impact parameter, in Fig. 4(b) we show 
with the star symbol the sensitivity parameter r evaluated by changing the option from 
Stiff 1 to Stiff2 and by only taking into account the contributions of the two main fragments. 
As we can see, the partial contribution shows a rather reduced sensitivity to the different 
options as compared to the case obtained by using the global information on the system. 

Finally, in the following we show the role of the correlation coefficient a, introduced in Sec. II, 
into determine the sensitivity of the investigated observable to the density dependence of 
the iso- vectorial interaction. For this aim, in Fig. 4(c), we compare, for different impact 
parameters, the values of r obtained for the Stiff2-Stiffl options with the ones obtained in 
the I.M.F.A. case. Fig. 4(c) clearly shows that in the I.M.F.A. case, at the investigated 
energies, the sensitivity of the isospin equilibration process to the behavior of the symmetry 
interaction is rather reduced. The I.M.F.A. also strongly affects the values of (V). In 
particular, independently from the used options, it produces for b r =0A, values of \(V Z )\ 
about four times larger than the ones obtained with full CoMD-II calculations indicating a 
reduced capacity to obtain isospin equilibration along the z direction. 

IV. SUMMARY AND CONCLUSIVE REMARKS 

In summary, in this work the isospin equilibration process has been investigated by study- 
ing the ensemble average of the time derivative of the total dipole (V) evaluated through 
CoMD-II calculations. Some general properties of this quantity have been discussed. In par- 
ticular, it allows to generalize the definition of isospin equilibration also in complex reactions 
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evolving through multi-fragmentation processes. As an example, calculations performed for 
the asymmetric charge/mass system 40 Cl+ 28 Si at 40 MeV/nucleon show that the asymptotic 
values of (v) for these processes are quite sensitive to different options for the iso- vectorial 
potential; moreover, in central and mid-peripheral collisions, the dipolar contribution asso- 
ciated to the pre-equilibrium emission of charged particles is relevant to determine the value 
of (v) and the related sensitivity to different density dependent form factors. Semi classical 
CoMD-II calculations performed in the so-called I.M.F.A. scheme also highlights the promi- 
nent role which could be played by the two-body neutron-proton correlations in the study 
of the dynamics leading to the isospin equilibration processes. 
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Fig. 1 - Average dipolar signals (T?) for 6=3 fm along the z direction are plotted as 
a function of time in the intervals 0-140 fm/c (panel (a)) and 110-300 fm/c (panel (b)). 
Different lines indicate different options for the iso- vectorial interaction (see the text). Panels 
(c) and (d) display the average dipolar signals along the x direction for 6=1.5 fm and in the 
same time intervals like panels (a) and (b) respectively. 

Fig. 2 - Asymptotic values (V x ) and (V z ) evaluated for different values of the b r parameter 
(see the text). Different symbols indicate different options for the iso- vectorial interaction. 

Fig. 3 - In panel (a) the Rql = |y§y ratio (see the text) is plotted as a function of the 
reduced impact parameter for the Stiff2 option. Panels (b) and (c) display the neutron- 
proton differential flow F np and the number of free nucleons Aq as a function of the cm. 
rapidity respectively. In the the panel (d), F np is plotted after the correction for the relative 
cm. velocity of the neutron and proton "gases" lf^ N . The vertical dotted lines represent 
the target and projectile rapidity. Different symbols refers to different options. In the panel 
(e) the average N ) contributing to the isospin equilibration process is represented for 
the Stiff 1 and Soft options. 

Fig. 4 -(a) (V x ) is plotted as a function of the corresponding (V z ) value for different 
reduced impact parameter b r values (b ma x — 7.5 fm and Ab r = 0.1 ) and different options 
(different symbols). The arrow indicates the direction of increasing impact parameters, -(b) 
relative changes r for the ratio R (see the text) evaluated for different couples of options as 
a function of b r . -(c) Values of r evaluated for the Stiff 1-Stiff2 options as a function of b r are 
plotted in the case of full CoMD-II calculations and in the case of I.M.F.A. approximation 
(see the text). The lines which join the points are meant only to guide the eye through the 
shown trend. 
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